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1. Introduction 

Primordial nucleosynthesis happens between 1 and 10^ seconds after the big bang at 
temperatures ranging from T ~ 1 McV to 70 keV. (These are the temperatures of weak 
freeze out and the end of the D bottleneck, respectively). Predictions of primordial light 
element abundances, D, '^He, "^He and ^Li, and their comparison with observations are a 
crucial test of the standard big bang cosmology. The uncertainties in these predictions are 
dominated by the nuclear physics input from reaction cross sections. Reaction databases 
are continuously updated [1, 2, 3], with more attention now paid to the error budget. 
In order to understand big bang nucleosynthesis (BBN) more clearly, it is essential to 
accurately measure these reaction cross sections at the energies relevant for BBN. 

The radiative neutron capture on a proton, np dj, is one of the key reactions for 
BBN, since this process is the starting point of the synthesis of the hght elements {i.e. it 
determines the end of the D bottleneck). The cross sections of the np reaction have 

been measured by Suzuki et al. [4] and Nagai et al. [5]. Its inverse process, the photo- 
disintegration of the deuteron, ^7 np, has had its cross section measured near threshold 
by Kara et al. [6] and Moreh et al. [7]^, and the photon analyzing power for the deuteron 
photo-disintegration are reported by Schreiber et al. [10] and Tornow et al. [11]. Although 
these data comprise an important data set, they nevertheless only sparsely sample the 
energies relevant for BBN. Hence an attempt at applying these experimental data directly 
to the BBN predictions would make the uncertainties larger. 

Theoretical calculations suggest that theory errors can be sufficiently smaller than 
typical experimental uncertainties ~ 5%, and that they can provide a very useful discrim- 
inant for theories and their perturbative schemes. BBN reaction compilations adopt these 
theory-based cross sections since they can provide more robust and accurate predictions 
than experiment alone. However, the uncertainties from the recent theoretical estima- 
tions of the cross section for np d'j at BBN energies are considerably different from 
each other; 4% [12], 2 ~ 3% [13], and 1% [14]. These differences could lead to different 
uncertainties in the BBN predictions, and thus it is necessary to examine the relevant 
error budget for the np d'-f process with a new method. 

Effective field theories (EFTs) provide a model-independent calculation and a system- 
atic perturbation scheme in terms oi Q/A in calculations of various low energy hadronic 
processes [15, 16], where Q is a typical momentum scale of a reaction in question and 
A is a large scale integrated out from effective Lagrangian. Since the energies relevant 
for BBN (T ~ 1 MeV) are significantly smaller than the pion mass, we can consider the 
pions as heavy degrees of freedom and integrate them out of the Lagrangian. Pionless 
EFTs [17] have been intensively studied in various two-, three- and four-nucleon processes 
for the last decade (see, e.g., Refs. [18, 19] for reviews and references therein). Conver- 
gence in the pionless EFT-based perturbative expansion turns out to be rather slow for 
the deuteron channel due to a relatively large expansion parameter, Q/A ~ 1/3 [20]. This 
large expansion parameter essentially determines the uncertainty estimates of the pion- 
less EFT calculations. For example, in the next-to-next-to-next-to leading order (N'^LO) 

^We do not include the data from Bishop et al. [8] in this work because of the wrong normahzation 
factor of the data [9]. 
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calculation of the np d'y cross section at the BBN energies, Chen and Savage found a 
(1/3)^ ~4% error [12]. Rupak pushed the calculation one order higher, i.e., up to N^LO, 
and found a (1/3)^ ~ 1% theoretical uncertainty in the cross section [14]. 

It has been suggested that the convergence of the pionless EFTs for the deuteron 
channel can be improved by adjusting the deuteron wave function so as to fit it to the 
asymptotic ^-state normalization constant Zd = jpii/il — ■jPd) (7 = V^nB where B is 
the binding energy of the deuteron and pd is its effective range) [21, 22, 23] . By introducing 
dibaryon fields which represent a resonance scattering state of the ^5*0 channel and the 
deuteron bound state of the ^5*1 channel, Beane and Savage showed that the dibaryon 
EFT (dEFT) without pions can naturally account for the long tail of the deuteron wave 
function in the renormalized dibaryon propagator at the deuteron pole [24].^ The slow 
convergence problem, however, was not fully resolved in dEFT. For instance, as discussed 
in detail in Ref. [28], when one includes the electromagnetic (EM) interaction into the 
np system in the framework of dEFT, vector(photon)-dibaryon-dibaryon (ydd) vertices, 
which are classified in the subleading order in the previous work [24], give contributions 
comparable to those of the leading ones. In Ref. [28], we suggested a simple prescription 
to extract a LO contribution from the low energy constants (LECs) of the Vdd vertices 
(we will mention it below), and this re-ordering of the Vdd term has shown a satisfactory 
convergence rate similar to that reported in other EFT calculations. We employ this 
modified counting of the Vdd vertex, and confirm that this machinery is useful for the 
calculations of processes and observables considered here. 

In this work, we calculate the cross sections of the np — > d^ process at BBN energies, its 
inverse reaction d'j — > np, and the photon analyzing power for the d'j — > np process within 
the framework of dEFT up to next-to leading order (NLO). Values and uncertainties of 
LECs that appear in the amplitudes are estimated by a Markov Chain Monte Carlo 
(MCMC) analysis using the relevant low energy experimental data; the total cross section 
of the np scattering at the energies <5 MeV, the rates of the np — > dj process, its 
inverse process d'j — > np, and the photon analyzing power in the d'j — > np process. 
Having fitted the values of the LECs, we compare our results with the experimental data 
mentioned above, and find that our statistical error bars of the np d'-f cross section 
are satisfactorily improved compared to the experimental ones. We also compare our 
result of the np — > d^ cross section with other theoretical estimations, the pionless EFT 
calculation up to N^LO [14], a calculation with Argonne vl8 (Avl8) potential and the 
meson exchange current [29], and the result of an R- matrix analysis [30]. Our result agrees 
quite well with those of the previous EFT and Avl8 potential model calculations within 
the uncertainties (~1%) estimated by MCMC. On the other hand we find significant 
(~4.6%) difference in the np d^ cross sections from the R-matrix theory estimated at 
£■ = 0.1 and 1 MeV where E is the total two-nucleon kinetic energy in the center of mass 
frame. Since the various theoretical uncertainties from higher order corrections in the 
pionless EFT calculation have already been studied in Ref. [14] and the slow convergence 

^Recently, dEFT has been employed in studying, e.g., neutron-neutron fusion [25], muon capture on 
the deuteron [26] , and the deeply- virtual Compton scattering dissociation of the deuteron and the EMC 
effect [27]. 
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problem in the former pionless EFT calculations has been resolved in dEFT, we discuss 
that the theoretical estimations of the np — > ^7 cross section at the BBN energies are 
reliable with an uncertainty of <1%. 

This paper is organized as follows. In Sect. 2, we present the pionless effective La- 
grangian with dibaryon fields up to NLO. We calculate the amplitudes with the S- and 
P-wave np states and the cross sections for the np — > (i7 and d'-) — > np processes and 
the photon analyzing power for the — > np process up to NLO in Sect. 3. Utilizing 
a Markov Chain algorithm, we determine the vahics and uncertainties of the LECs in 
Sect. 4. In Sect. 5, we compare the resulting obscrvables to the experimental data and 
other theoretical estimations. In Sect. 6, we summarize our results and give discussion. 
In appendix A, we present the expressions of the renormalized dibaryon propagator and 
the S-wme NN scattering amphtudes, and show the relations between the LECs in the 
strong sector and the parameters in the effective range theory. In appendix B, we describe 
in detail the MCMC analysis in determining the LECs and cross sections. 

2. Pionless effective Lagrangian with dibaryon fields 

A pionless effective Lagrangian for micleon and dibaryon fields interacting with an 
external vector field can be written as [24, 28] 



st ■) 



where is the nucleon Lagrangian, Cg and Ct are the Lagrangian for the dibaryon fields 
in ^5*0 and ^Si channels, respectively. Cgt is the Lagrangian that accounts for the isovector 
EM interaction of the dibaryon fields. 
>Cjv in the heavy-baryon formalism reads 



jC^^N^[zvD + ^ [{v ■ Df -D'- i[S^, S'^] ((1 + Ky)f+ + (1 + Ks)f,^)] I AT , (2) 



where is the velocity vector satisfying v"^ = 1; we take v'^ = (1,0). iS^ is the spin 
operator 25*^ = (0,ct). = — ■ — fV^, where and are the external 

isovector and isoscalar vector currents, respectively. = | • {df^Vu — duVn) and f^^, = 
^{d/jV^ —d^V^). niN is the nucleon mass and (ns) is the isovector (isoscalar) anomalous 
magnetic moment of the nucleon, Ky — 3.706 {ks — —0.120). 
jCs, C^ti and Lst for the dibaryon and two nucleon fields read 
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The covariant derivative for the dibaryon field is given by = — iCV^^^, where V^^* 
is the external vector field and C is the charge operator of the dibaryon fields; C = 0, 1, 2 

— * 

for nn, np, pp channels, respectively, and wc have set e = 1. B is the magnetic field given 
by i? = V X V^^*. at (o"s) is the sign factor, (A^) is the difi^erence between the dibaryon 
mass rrit (m^) in the ^5*1 (^5*0) channel and the two-miclcon mass; rrit^s = 2mjv + A^^^, 
and Us^t are the dibaryon-nucleon-nucleon (dNN) coupling constants of the dibaryon spin 
singlet and triplet channels. In appendix A, these LECs in the strong sector are related 
to the parameters of the effective range theory, li and I2 are LECs of the Vdd vertices 
and can be fixed, for instance, by the thermal np d'-f rate and the deuteron magnetic 
moment, respectively. Wc note that wc have separated the leading contributions from 
the coefficients in the Vdd vertices associated with li and I2 and fixed them mainly 
by one-body interactions, i.e., the vector (photon)-nucleon-nucleon (VNN) couphngs, as 
suggested in Ref. [28]. Consequently, h and I2 terms in this work give genuine NLO 
contributions, pd and tq arc the effective ranges for the deuteron and ^Sq scattering state, 
respectively. P/'^'' is the projection operator for the S = '^'^^^Lj channel. For the S and 
P waves which are dominant at low energies, the projection operators are given as [31]; 

= ff'P'-^a^r^r^ > / ^Tr(pf )t^(5)) ^ 1^^^.^ 

where e* and e*-' are J=l and 2 polarization tensor, respectively, and ai (tq) with i{a) — 
1,2,3 is the spin (isospin) operator. 

We adopt the standard counting rules of dEFT [24]. Introducing an expansion scale 
Q < A (~ m7r), we count magnitude of spatial part of the external and loop momenta, |p| 
and |/[, as Q, and the time component of them, p° and P, as Q^. Thus the nucleon and 
dibaryon propagators are of and a loop gives a factor of due to the 4-dimensional 
differential volume in the loop integration. The scattering lengths and effective ranges are 
counted as Q ~ {7, I/oq, 1/pd, l/''^o}- This ensures, as discussed in the introduction, that 
one reproduces the long tail of the deuteron wavefunction characterized by -\/^ thus 
has good convergence [21, 22, 23]. Orders of vertices and diagrams are easily obtained by 
counting the numbers of these factors. 

3. Amplitudes and cross sections 

Diagrams for np d'-f up to NLO are depicted in Fig. 1. Diagrams (a) and (b) give 
only LO contributions. Diagram (c) consists of LO and NLO fractions and the NLO 
contributions from the diagram (c) stem from the LECs li and I2, as discussed above. 

Summing up the contributions of the diagrams (a), (b) and (c) in Fig. 1, we obtain 
the amphtudes for the initial ^Si and ^Sq states as [24, 28] 

iA^^+'^+'^\'S,) = -^(e^,) X el) . (e7^) x k) 
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(a) (b) (c) 

Figure 1: Diagrams contributing to np —>■ and —>■ np: Diagrams (a) and (b) are 
of LO, 0{Q^/'^), while diagram (c) consists of LO and NLO terms, ^(Q^/^). The single 
lines, the double lines with a filled circle (see Fig. 6 as well), and the wavy lines denote 
nucleons, renormalized dibaryons, and photons, respectively. VNN vertex "x" in (a) and 
(b) and the LO part of Vdd vertex "x" in (c) are proportional to (1 + k,s) and (1 + /ty) 
for the initial ^5"! and ^So channel, respectively. LECs h and I2 appear in the NLO part 
of the Vdd vertex "x" in the diagram (c). 
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Here 2p is the relative momentum of the two-nucleon system, and k is the momentum of 
the outgoing photon; p = \p\, k = \k\, and k = k/k. e^^^-j and e^^-, are the polarization 
vectors for the outgoing deuteron and photon, respectively, is the scattering length in 
the ^Sq channel. The LECs in the strong sector have been renormahzed by the effective 
range parameters. See appendix A for details. We note that the isovector Ml amplitude 
(the ^5*0 channel), Eq. (8), has contributions from both LO and NLO. Whereas, the 
isoscalar Ml amplitude (the ^Si channel), Eq. (7), has no LO contribution due to the 
orthogonality between the bound and scattering states for the ^5*1 channel, but non-zero 
amplitude proportional to I2 appears at NLO. 

From the diagram (a) in Fig. 1, amplitudes with initial P- waves are obtained as; 



I. 



u(«)('A) = -i6i • (67,) X 7c^))\/t^\/I r L 2w (10) 

•yl(a)/3p>, _ ijM *j / ^^7 2 p 

Ki) w\j I- ^3 mjv(7 + P ) 
We note that the diagrams (b) and (c) in Fig. 1 do not contribute to the P-wave amplitudes 
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because the states of the renormahzed dibaryon propagator are only the ^-waves. 

Having the amphtudes calculated above, we can easily obtain the expressions of the 
cross sections of the np — > d'-y and ^7 — ^ np processes and the analyzing power for the 
(i7 — s> np process. The total cross section of the np — > d'y process in the CM frame 
reads [24, 28] 

^^a(Tl±P^^I^|,_ (12) 

" spin 

where a is the fine structure constant. ^ 

The total cross section of the photo-disintegration of the deuteron, 1^7 np, has a 
simple relation with the cross section of its inverse process as [12] 

(7(7^ np) = ^\ {np d-f) , (13) 

where is the photon energy. Since we have already obtained the np-capture cross 
section, the calculation of a^'-fd np) is straightforward. 

The photon analyzing power S(6') with the linearly polarized photons in the d^ — > np 
process is defined as E(6') = {N\\ — N±)/{N\\ + N±) where N\\ and A^_l are the number 
of outgoing neutrons counted in and out of the horizontal 7-ray polarization plane, re- 
spectively, and 6 is the angle between the incoming photon and outgoing neutron in the 
laboratory frame. This quantity is related to the Ml and El contributions to the total 
cross sections of the dj — > np process, ctmi and aEi, respectively. The relation can be 
found e.g., in Eq. (4) of Ref. [10], which reads 

m - (14) 

gmi and gei are easily calculated by using Eqs. (12, 13) and the expressions of the Ml 
and El amplitudes obtained in Eqs. (7, 8, 9, 10, 11). 

4. Parameter determination from experimental data 

In this section, we determine the values and uncertainties of the LECs that appear 
in our results from the relevant low energy experimental data. The cross sections and 
photon analyzing power obtained from the amplitudes in Eqs. (7, 8, 9, 10, 11) depend 
on six parameters; Oq, Tq, 7, Pd, /i, and 12- 7 and I2 are precisely determined by the 
deuteron binding energy B and magnetic moment ^a-i respectively. Using the relation, 
jid — ^ -\- i^s -\- Zdk: we have I2 — —0.0154 fm. It is ~50 times smaller than the value of 

''In obtaining the total cross section, the following identities are useful. 

spin spin spin spin 

We ignore the amplitude for the channel in Eq. (7), as will be discussed later. 
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li in Table 1 (and the expressions for tlic li and I2 terms in the amphtudes are almost 
identical). We will consider the error of the cross section at the order of 0.1%, so we 
neglect the amplitude in Eq. (7) from the ^5*1 channel in the following calculations. Thus, 
the parameters a^, Tq, p^, and li are determined from the low energy experimental data, 
whereas 7 is fixed from B and treated as a constraint in the fitting. 

We fit these parameters to the low energy np data with a Markov Chain Monte Carlo 
(MCMC) analysis using the Metropohs algorithm [32]. This method is more efficient in 
exploring parameter space than creating a mult i- dimensional grid of parameter values 
and interpolating to find the underlying likelihood distribution and favored parameter 
values. The method is a random walk constrained to "walk" in regions of low x^- As 
the increases for a particular step, the probability of the step being accepted into the 
chain decreases. This "walk" explores the allowed parameter space, and if long enough 
will explore all channels of parameter degeneracy. The resulting sample, then provides a 
direct probe of the parameter likelihood, and can be used to determine, e.g., the means, 
standard deviations and correlations of the parameters. It also provides an accurate way 
of propagating the uncertainties and correlations of the parameters into the cross section 
uncertainties, which depend non-linearly on the parameters. One drawback is that this 
method is not particularly efficient at sampling multi-peaked distributions. The random 
walk can become trapped in a local minimum. For the case we are considering, the 
5-dimensional likelihood distribution is singly peaked {i.e. a unique solution exists). 

We now list the low energy experimental data employed in our fitting. To include 
the high precision measurements of the deuteron binding energy, we adopt the value of 
B = 2.2245671 ±0.0000042 MeV^ [34] as the additional constraint mentioned above. This 
accurate value of B gives the uncertainty of 7 in the order of 10~^, thus it is effectively 
fixed, independent from the other constraints. The np-scattering data are found on the 
NN-Online website [35]. We restrict ourselves to the angle- integrated scattering cross 
sections with center-of-mass energies <5 MeV, and have Ngct — 2124 data points. We 
adopt the 2 thermal neutron capture cross sections from Cox et al. [36] (334.2±0.5 mb) 
and Cokinos and Melkonian [37] (332.6 ± 0.7 mb),^ and also adopt the np-capture data 
of Suzuki et al. [4] and Nagai et al. [5] , the photo-dissociation cross section data by Hara 
et al. [6] and Moreh et al. [7], and the photon analyzing power data from Schreiber et 
al. [10] and Tornow et al. [11]. 

Using the experimental data mentioned above, we fit these parameters employing the 
MCMC analysis. The steps in doing the MCMC analysis are described in detail in ap- 
pendix B. We initialize the Markov chain at the point in 5 dimensional parameter space 
that minimizes the . After verifying that the had only one minimum, we did not 

^The quoted deuteron binding energy is a weighted average of the available measurements. The 
error is the weighted dispersion about the mean, recommended by [3, 33] because standard techniques 
underestimate uncertainties when data arc discrepant. 

^In a private communication mentioned in [38], the Cox et al.'s thermal cross section is renormalized 
from 334.2 to 333.9 mb. This likely indicates hidden systematics in the experiment and subsequently 
larger (and unacknowledged) uncertainties, though we adopt the original published number. In addition, 
the thermal np capture rate, a = (334 ±3) mb, has been estimated by T.-S. Park et al. in the calculation 
of pionful effective field theory [39]. 
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Figure 2: Contour plots of two dimensional parameter spaces for 5 dimensional probability 
distribution generated by the MCMC analysis. We draw the contours of 1,2, and 3 a for 
each set of the parameters, ao, ro, Pd,, and h, whereas B is effectively fixed at the observed 
value and uncorrelated with the other parameters, and thus not shown. 



need to run more than one Markov chain, making the analysis much more efficient. We 
display, in Fig. 2, two dimensional contour plots of 1, 2, and 3 a for each set of param- 
eters, ao, ro, Pd, h from our 5 dimensional probability distribution generated from the 
MCMC analysis. We note that since B is effectively fixed and uncorrelated with the other 
parameters, contour plots with B are not necessary. "MCMC" in Table 1 are the average 
parameter values and standard deviations from this analysis with the total number of the 
data A^tot = 2147 and the minimum x^, Xtot miu ~ 2303.00. (For more details, see Table 4 
in the appendix B.) These parameter values agree well with prior determinations shown 
in "Prev. Meth." in the same table. The small (< 2%) differences between "MCMC" 
and "Prev. Meth." could be interpreted as effects from either higher order corrections not 
included in our calculation or simply statistical fluctuations. Also in Table 1, we show the 
correlations of the parameters. We find weak correlations of with the three parameters 
ro, Pd, and h and strong correlations among these three parameters. Though we formally 
counted the four parameters as in the same order in this work, this may indicate the 
existence of a perturbative series that the contribution of ao is the LO one and the three 
parameters are in the same order (NLO), as already discussed in Ref. [20]. 

5. Numerical results 

First, we present our numerical results obtained by using the fitted values of the pa- 
rameters "MCMC" in Table 1, and compare them with the experimental data relevant to 
the BBN energies. We plot in Fig. 3 the total cross section (in mb) of np d'-f multiplied 
by the speed (in m/ns) of the neutron in the laboratory frame as a function of the inci- 
dent energy of the neutron En- We also plot the Ml and El contributions to the total 
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1.000 
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1.000 


0.898 
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0.782 ± 0.022 


0.893±0.038 


0.496 


0.935 


0.898 


1.000 



Table 1: Values and correlations of the parameters. The values of the parameters obtained 
are in units of fm. "Prev. Meth." are the adopted parameter values from previous 
works; of these ao, ro, and p,j arc from Rcf. [40] and li from the thermal np capture 
rates [28]. "MCMC" is obtained by the MCMC analysis (x^-fit) by using the low energy 
experimental data. Also shown are the parameter correlations. They are defined by 
COR{i,j) = COV"(i, j)/((Tstat.(0'^stat. (i)), whcrc COV{i,j) is the covariance between the 
i^^ and j*^ parameters and (Jstat. (i) is the statistical uncertainty of the i^^ parameter. 



capture cross section in Fig. 3, where the Ml contribution comes from the amplitude of 
the initial ^5*0 state in Eq. (8) and the El contribution from the amplitudes of the initial 
P-wave states in Eqs. (9,10,11). At very small energies, the Ml contribution overwhelms 
the El one. They become similar at around En ~ 0.45 MeV and after that, the cross 
section is dominated by the El contribution. ^° The experimental data of np — > d'j are 
compared with our result in Fig. 3. Our results lie within the errors (5~6 %) of all the 
data by Suzuki et al. and Nagai et ai. In Fig. 4, we plot our result of the cross section of 
the dj — > np process and also separate the contributions from the Ml and El transition 
amphtudes. Recent measurement of the cross section at the BBN energies is reported by 
Kara et al. [6] . An old datum by Moreh et al. [7] is also included in the figure. Our result 
agrees well with these experimental data within the error bars. In Fig. 5, we plot our 
results of S(^) at = 90° and 150° where experimental data are available [10, 11]. We 
find good agreement between our predictions and the experimental data. The error bars 
estimated for these quantities, discussed in detail below, are much improved compared to 
those of the experimental data. 

Now, we compare our results with the predictions of the np d'-f cross sections at 
the BBN energies from various theories in Table 2. Our results, "MCMC" and "Prev. 
Meth.", in Table 2 are calculated from the amplitudes in dEFT up to NLO by a MCMC 
analysis and by using the values of the "Prev. Meth." parameters in Table 1, respectively. 

^°Tlie contributions of Ml and El transitions to the total cross section of np dj and d^y np are not 
disentangled in the measurements, but the role of each amplitude can be studied from the measurement 
of 2(6*). 

^^One may notice a departure (more than 1%) of our estimation from the center values of experimental 
data in Figs. 3 and 4. This is because the curves plotted in the figures are mainly determined by the 
other data: e.g., the np scattering data have a prime role to determine the energy dependence of the 
curves, whereas the accurate thermal np capture rates determine the normalization of them. However, 
as verified by the good x^, the curves are well within the experimental error bars. 
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Figure 3: Total cross section ofn + p— >(i + 7in units of mb multiplied by neutron speed 
in m/ns. Neutron energy is in the laboratory frame. Dashed and dotted curves are 
the Ml and El contributions to the total cross section, respectively. The experimental 
data are from Suzuki et al. [4] and Nagai et al. [5]. 



E (MeV) 




Cross section 


^mb) 




MCMC 


Prev. Meth. 


Rupak[14 


Nakamura[29 


Hale [30 


1.265 X 10-*^ 


333.8(15) 


333.7(15) 


334.2(0) 


335.0 


332.6(7) 


5 X 10-^ 


1.667(8) 


1.666(8) 


1.668(0) 


1.674 


1.661(7) 


1 X 10-^ 


1.171(5) 


1.171(5) 


1.172(0) 


1.176 


1.167(2) 


5 X IQ-^ 


0.4979(23) 


0.4976(21) 


0.4982(0) 


0.4999 


0.4953(11) 


1 X 10^2 


0.3322(15) 


0.3319(14) 


0.3324(0) 


0.3335 


0.3298(9) 


5 X 10-2 


0.1079(5) 


0.1079(4) 


0.1081(0) 


0.1084 


0.1052(9) 


0.100 


0.0634(3) 


0.0634(2) 


0.06352(5) 


0.06366 


0.0605(10) 


0.500 


0.0341(2) 


0.0343(1) 


0.0341(2) 


0.03416 


0.0338(8) 


1.00 


0.0349(3) 


0.0352(1) 


0.0349(3) 


0.03495 


0.0365(8) 



Table 2: Theoretical predictions of the total cross section of the n + p — > d + 7 process at 
the BBN energies. E is the energy of two nucleons in the center of mass frame. See the 
text for details. 
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Figure 4: Total cross section for the + 7 — n + p process. Dashed and dotted curves 
represent the Ml and El contributions to the total cross section, respectively. The 
experimental data are from Hara et al. [6] and Moreh et al. [7]. 
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Figure 5: Photon analyzing power S(6') for the (i + 7 — n + p process. We plot our results 
at ^ = 90° and 150°. The experimental data are from Schreiber et al. {9 — 150°) [10] and 
Tornow et al. = 90°) [11]. 
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Wc note that the values of the cross section in the column "MCMC" are the most likely 
cross section values with 68.3% central confidence limits. Thus we find that the error bars 
for the cross section are < 1%^^. Values obtained by Rupak [14] in Table 2 are from the 
pionless EFT (without dibaryon fields) up to N^LO and those by Nakamura [29] are from 
the potential model calculation using wave functions from the Argonne vl8 potential and 
meson exchange currents. The results by Hale are obtained from an R- matrix analysis [30] . 
There is good agreement between our two results, "MCMC" and "Prev. Meth." , up to 0.1 
MeV, whereas they show small differences (■-^0.6% and 0.9%) at E = 0.5 and 1 MeV. We 
also find good agreement between the "MCMC" analysis with that of the pionless EFT 
up to N^LO by Rupak (< 0.2%) and with that of the accurate potential model including 
the exchange current by Nakamura (< 0.5%), while the results of the R-matrix theory at 
E — 0.1 and 1 MeV significantly differ from the other estimations by ~4.6%. 

Finally, we determine and present a thermal capture rate and relative error taking the 
recommended cross section from "MCMC" in Table 2 for use in BBN computer codes: 

/ = NA{av) = 44216.0[cm^ s"^ g-^](l + 3.75191 Tg + 1.92934 Tg^ 
+0.746503 r| + 0.01970237;^ + 3.00491 x 10"^ r|) 
/(I + 5.4678 Tg + 5.62395 Tg + 0.489312 r| + 0.00747806 T^) , (15) 
5f/f = 0.00449213(1 + 3.08947 Tg + 0.13277 Tg^ + 1.66472 r|) 

/(I + 2.75245 Tg + 1.40958 Tg^ + 0.8791 T|) , (16) 

where Tg is a dimensionless temperature defined by Tg = T/{10^K). Using the thermal 
rates for the np-capture reaction from "MCMC dEFT", Rupak, and Hale, we show in 
Table 3 how the light element abundance predictions vary with the different np-capture 
cross sections. We find little change for the mass fraction of ^He, Yp and the mole fractions 
'^Hc/H, a tiny (~ 0.6%) change for D/H, and a small (~ 2.9%) change for ^Li/H in the 
predicted light element abundances. Differences between the light element predictions are 
not significant compared to the current estimated errors in the BBN predictions. In fact, 
the error budget in the BBN predictions is dominated by the errors in other reactions 
such as d{p, 7)^He and ^He(Q;, 7)^Be (see, e.g., Ref. [3] for more details). We have verified 
that the np-capture rate is not yet a major source of uncertainty in the light element 
abundance predictions. 

6. Conclusions and discussion 

In this work, we calculated the total cross sections for the np d'-f and d'-f np 
processes and the photon analyzing power in the dj — > np process at the energies relevant 
to BBN. The pionless EFT that incorporates dibaryon fields was employed, and the tran- 
sition amplitudes were calculated up to NLO. The values of the parameters (equivalently 
the LECs in effective Lagrangian) and their uncertainties are evaluated by a MCMC anal- 
ysis (x^-fit) using the relevant low energy experimental data. Comparing our results with 
the experimental data and the previous theoretical estimations, we find good agreement 

^^One should note that the error bars of the 68.3% confidence hmits are different from (and larger 
than) those of the standard deviation in multi-dimensional fits, whose typical error bars are < 0.3% [41]. 
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np-capture rate 




D/H X 10^ 


=^He/H X 10^ 


^Li/H X 10^" 


MCMC dEFT 


0.24852 


2.5467 


10.0921 


4.4646 


Rupak [14] 


0.24853 


2.5434 


10.0920 


4.4902 


Hale [13] 


0.24849 


2.5580 


10.0864 


4.3632 



Table 3: The table shows how the light element abundance predictions vary with different 
np-capture cross sections. The nuclear rate compilation from [3] was adopted, varying 
only the np-capture rate. In fact, these predictions follow the simple abundance scalings 
(Eqs. (44-47) in Ref. [3]), using relative values of the np-capture cross sections at the end 
of the D bottleneck {E ~ 0.07 MeV). 



within the estimated uncertainties (<1%) except for the np-capture rate estimated by the 
R-matrix analysis at £^ = 0.1 and 1 MeV where E is the initial energy for two-nucleon 
in the CM frame. These two values estimated by the R-matrix theory are considerably 
different from the other theoretical estimations by ~ 4.6%. Therefore, it would be impor- 
tant to experimentally measure the np-capture cross sections at these energies to resolve 
this significant discrepancy. 

Now we are in the position to discuss the uncertainties of np — > d-^ cross section at 
BBN energies. As discussed earlier, the EFT calculations provide model-independent 
expression for the amplitudes with a systematic perturbative scheme. In the pionless 
EFT calculation up to N^LO by Rupak, various corrections in the higher order terms 
have been taken into account and it has been concluded that the theoretical uncertainty 
up to the N^'LO calculation in the np — > d-^ cross section is less than 1%. The effective 
range corrections are resummed in dEFT, and so the convergence of dEFT is better than 
the pionless EFTs. Thus we expect that the contributions from the higher terms, i.e., 
the theoretical uncertainties in our calculation, to be less than Rupak's estimation. With 
the overall good agreement between our MCMC analysis and the Rupak and Nakamura 
calculations, the conclusion that the theoretical uncertainties in the np — > d-^ cross sections 
is <1% is well justified. The disagreement with the R-matrix analysis could be inferred 
from problems that R-matrix theory has in describing non-resonant reactions. 

Most of the recent theoretical calculations, as seen in Table 2, predict similar results 
with similar uncertainties, and the accuracy of the calculations is better than the np- 
capture experimental results at present. This is the case mainly because there are a lot of 
accurate experimental data for the np scattering and, we could accurately constrain the 
four effective range parameters from them. We had only one additional LEC li to fit from 
the low energy np d'-f and dr/ np cross sections and the photon analyzing power for 
the d'y np process. This situation, however, will change once we start studying other 
processes involving more than two nucleons for BBN, facing significant model dependence 
and lack of experimental data. As discussed in Ref. [42], the EFT approaches would be 
useful in studying few-body nuclear astrophysical processes because it provides simple 
model-independent expressions of the amplitudes with a finite number of LECs as well 
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Figure 6: Diagrams for "dressed" dibaryon propagator at leading order. A double (single) 
line stands for a dibaryon (nucleon) field. 

as a systematic expansion scheme. ^'^ We expect that the combination of dEFT and the 
MCMC analysis can be a useful tool to estimate reliable uncertainties of few-body nuclear 
reactions for BBN with the aid of relevant low energy experimental data. 
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Appendix A 

In this appendix we review the derivation of the renormalized dibaryon propagator 
from the low energy 5*- wave NN scattering [24, 28]. 

The leg's as,t and ys,t in the effective Lagrangian, Eqs. (3) and (4), can be fixed from 
the effective range parameters of the np scattering in ^Sq and states. Firstly, we derive 
"dressed" dibaryon propagators. LO diagrams for the dressed dibaryon propagators in 
the S wave channels are depicted in Fig. 6. Since the insertion of the two-nucleon one- 
loop diagram does not alter the order of the diagram, the two-nucleon bubbles in the 
propagators should be summed up to infinite order. Thus the inverse dressed dibaryon 
propagators for the spin singlet (s) C'So) and triplet (t) {^Si) channels in the center-of- 
mass (CM) frame read 



where we have calculated the two-nucleon one-loop diagram using the dimensional regu- 
larization. p is the magnitude of the nucleon momentum in the CM frame, and E is the 
total energy E ~ p^/mAr. 

^^Recently, the nd process has been studied in the pionless dEFT [43]. 
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Figure 7: Diagram for the NN scattering amplitudes. The dNN vertex is proportional 
to ys,t and the propagator of the dressed dibaryon field (a double line with a filled circle) 
is obtained from the diagram in Fig. 6. 



The S-wave NN scattering amplitudes for both spin channels obtained from Fig. 7 
read 

4:7r i 

iAs,t = {-iys,t) {iD,^tip)){-iys,t) = i^^TTATl S^F^T^^ , (18) 

and they are related to the ^'-matrix via 



Sst-l^ e'"'-'* - 1 = ^ = i ^--^ Mat , (19) 

where 5s,t are the 5'-wave phase shifts. The effective range expansion reads 

11 1 

pcoi5s = h -rop^ H , pcot6t = -7 + -pd(p^ + 7^) H , (20) 

do z z 

for the ^Sq and ^Si channel, respectively. Comparing the expressions of the amplitudes 
in Eqs. (18) and (19), one has as,t = — 1 and 

ys = — J— , Ds{p)^— — r— 1 — ^, (21) 

rriN^ pd 2 -f + ip- ^pd{p^ + 7 ) E + B 

where Zd is the wavefunction normalization factor of the dcuteron around deuteron bind- 
ing energy B. Ellipsis denotes corrections that are finite or vanish at E = —B. Thus one 
has 

7Pd ^23) 



1 -7Pd 

which is the same as the asymptotic 5*- wave normalization constant. 

Appendix B: Running a MCMC 

In this appendix we describe the steps in doing the MCMC analysis for determining 
the D = 5 parameters and cross sections. 
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1. Localizing the chain by minimizing the x^. 

The evolution of a Markov chain is governed by the values at various points in 
parameter space. We adopt the standard definition: 



where (Ti(expt) and (5crj(expt) are the experimentally measured values (^e.g. cross sections) 
and their total errors, respectively, and (jj(thry) is the evaluated theoretical value; the 
sum is over all data points. 

Before we evaluate the x^, we should pick some reasonable model space. Generally 
relying on previous works or positive definiteness, one can place limits on the allowed 
parameter space. Of course one can expand this model space if a minimum is found at 
an edge of the parameter space. 

We begin by selecting an initial set of parameters (po) and estimate step sizes (5po) 
be some fraction of the size of each parameter space direction. Drawing a set of random 
numbers z*with zero mean and unit variance {e.g. a Gaussian-normal random number), 
one can determine a new set of parameters via the relation: 



Evaluating the x^(p) of this new set of parameters, and comparing to the original X^(Po) 
we can determine if the new parameter values better describe the data. If the new x^ 
is smaller we accept the new parameter values (po — V)-, otherwise keeping the original 
parameter values. We calculate new parameter values and repeat. The parameter values 
will gradually evolve to the minimum possible x^i by incrementally decreasing the step 
size one can determine the best fit to some desired accuracy. 

2. Checking for convergence to unique minimum. 

By repeating this procedure with different starting parameter values in our model 
space we can determine the uniqueness of this minimum. This is particularly useful in 
testing the boundaries of the chosen model space. If a minimum is found on a boundary 
the parameter space must be enlarged. If any 2 parameters are completely correlated, a 
unique minimum will not be found {e.g. h and I'^y^ and one must reconsider the allowed 
parameter space. We will assume from now on that there is only one minimum in our 
model space and that no two parameters are completely correlated. 

3. Determining an appropriate step-size. 

To make the MCMC as efficient as possible, one needs to determine an appropriate step 
size, 6p. A simple method of estimating this is by varying individual parameter values 
away from the minimum, until the difference x^ ~ Xmin is unity. This choice makes for 
a good first estimate, and will generally be smaller than or equal to the true parameter 
errors. 

4. Running the chain(s). Enforced constraints. 

^^I'l is a LEG associated with a vector-dibaryon-nucleon-nucleon vertex. The LECs li and Z; are almost 
redundant in the np-capture cross section [28]. 




(24) 



p^po + 5poz . 



(25) 
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There are many ways to initiate a MCMC Some choose to pick an initial point ran- 
domly in the model space, while others choose to initiate the chain at the minimum. The 
latter method avoids the dependence on the prior probability distribution the initial point 
is generated from and thus reducing the overall convergence length of the chain. 

Once the initial point of the chain is chosen, we follow a procedure quite similar to 
that used in the minimization algorithm. One generates a new set of parameters just as 
in the minimization algorithm. There is then a set of criteria for accepting or rejecting 
the new parameter set: 

1. if Ax^ — X^iP) ~ X^(Po) < we accept the new parameter values. 

2. otherwise there is a finite probability, P = exp (— for accepting the point. 

For a reasonable step size choice, this allows for the efficient exploration of the "tail" of 
the parameter likelihood distribution. Whether or not one accepts the new point, a new 
set of parameters is drawn and this step is repeated until the chain has met its convergence 
criteria or some maximum length. 

5. Checking convergence to a "full" sample. 

A relatively simple method to check the convergence of a single chain is to calculate 
the 1** and 2"*^ order moments of the chain of a specific length N. We thus calculate the 
means, variances and covariances of the D parameters: 

1 ^ 

m = -^Pi, (26) 

i 

c{N) = -J;^J:\Pi-m)]^\p^-m, (27) 

i 

where the variances are the diagonal components of the covariance matrix. As the MCMC 
converges, these moments of the underlying likelihood distribution will plateau and the 
fractional difference between the A^*'* and (A'+l)*'^ moments (or functions of the moments, 
e.g. det(C)) should fall like ~ 1/A". Thus, we choose N in such a way to reach a certain 
desired fractional uncertainty. 

One can also compare multiple chains and their moments. One can then compare the 
variance of a single chain to the variance of the chain means, adopting the convergence 
criteria from [44]. However, since we are starting our chains at the minimum, a single 
chain is all that is needed once the chain has grown longer than the intrinsic correlation 
length of the chain (typically ~ few 100 's steps) and the convergence criteria for the single 
chain is all that is needed. 

To meet convergence criteria, chains with length N ~ 10^ are required. The parameter 
step size is adjusted so that the acceptance is between 25 and 50%. We find that a 40% 
acceptance rate is most efficient. 
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